On the usage of Flaring Gas Layers to determine the 
Shape of Dark Matter Halos 

To appear in the August 1995 issue of The Astronomical Journal 



Rob P. Oiling 

Columbia University 



ABSTRACT 

I present a new method of deriving the shape of the dark matter (DM) halos of spiral 
galaxies. The method relies on the comparison of model predictions with high spectral 
and spatial resolution HI observations of the gas layer. So far, determinations of the 
flaring of the gas layer (i.e. the increase of the thickness with galactocentric radius) 
have been used to determine the mass-to-light ratio, Ai/C , of the stellar disk of several 
edge-on galaxies. In this paper I describe a method which can be used to determine 
the shape of DM-halos. This technique will be applied in a forthcoming paper. 

I show that the model predictions of the gas layer width are best calculated using 
a global approach, in which the potential arising from the total mass distribution of 
the galaxy is used in the calculation of the vertical distribution of the gas. I developed 
a new algorithm to calculate the force field of an arbitrary, azimuthally symmetric, 
density distribution. This algorithm is used to calculate the forces due to the radially 
truncated stellar disk as well as of the fiaring gas layer. 

I use a simple two-parameter family of disk-halo models which have essentially the 
same observed equatorial rotation curve but different vertical forces. This mass model 
is composed of a stellar disk with constant , and a DM-halo with a given axial 

ratio (Sackett & Sparke [ApJ, 1990, 361, 408]). I approximate the radial force due to the 
gaseous disk, and iteratively determine the vertical force due to the global distribution 
of the gas. 

In agreement with Maloney [ApJ, 414, 41, 1993] I find that beyond the Holmberg 
radius, the thickness of the gaseous disk is sensitive to both the fiattening of the DM- 
halo and the self-gravity of the gas. I also show that the inferred DM-halo fiattening is 
not sensitive to the particular choice of disk-halo decomposition. 

I show that the determination of the thickness of the gas layer is not restricted to 
edge-on galaxies, but can be measured for moderately inclined systems as well. Thus, 
in combination with detailed modeling, high resolution HI imaging of nearby galaxies 
with extended HI envelopes will enable us to determine the shape of the DM halo of 
these galaxies. 
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1. Introduction 

Until recently, little was known about the shape of 
dark halos : measurements of the equatorial rotation 
curve provide a one dimensional probe to the poten- 
tial only. It is possible to construct both spherical 
and flat mass models which generate the same rota- 
tion curve. Therefore, the shape of the dark halo can 
not be inferred from galactic rotation curves alone. 

Existing methods to estimate the shape of DM- 
halos have given mixed results. It has been suggested 
(Dekel & Shlosman, 1983; Toomre, 1983; Sparke & 
Casertano, 1988) that galactic warps, frequently ob- 
served in the outer parts of gaseous disks, result if 
the galactic disk is tilted with respect to the plane of 
a flattened dark halo. Hoffner & Sparke (1994) in- 
vestigated the time evolution of such warps and con- 
cluded that only one out of the five systems they 
studied requires a halo as flattened as E6. Polar 
ring galaxies, early type spirals with material rings 
in a plane perpendicular to the galaxy plane, probe 
the dark halo potential in two perpendicular direc- 
tions. Careful analysis of the stellar and gas dy- 
namics of the polar ring system NGC 4650A, led 
Sackett & Sparke (1990, hereafter SS90) and Sackett 
et al. (1994, SRJF94 hereafter) to conclude that the 
dark halo in this system could be as flattened as E6- 
E7 ((?=0.4 to 0.3). Recently, Sackett et al. (1994b) 
reported the discovery of a rather flattened (q w 0.5) 
luminous halo around the edge-on spiral NGC 5907, 
which can account for the observed rotation curve 
provided that this luminous material has a large mass- 
to-light ratio (w 450). 

Like the polar rotation curve, the thickness of the 
gas layer depends on the vertical force, Kz^ and hence 
on the shape of the DM-halo (Maloney 1992; Olhng & 
van Gorkom 1992; Kundic et al. 1993, KHG93; Mal- 
oney 1993, M93). That the shape of the dark halo 
influences the vertical distribution of the gas can be 
easily understood. Consider a round halo, with a cer- 
tain density distribution, then decrease the vertical 
scale-height while maintaining the total mass. Con- 
sequently, the densities as well as the exerted grav- 
itational forces will increase, resulting in a thinner 
HI disk and higher rotation speeds. In order to keep 
the same equatorial rotation curve, one has to de- 
form the DM-halo in a very specific way (Appendix 
A, Figure ||) as a result of which the DM-halo densi- 
ties (at large distances) will be roughly inversely pro- 
portional to the flattening q (= c/a). In accordance 



with Maloney (1993), I find that (at large distances) 
the thickness of the gas layer is roughly proportional 
to the square root of the halo flattening, q. Several 
authors have employed measurements of the gas layer 
thickness to infer the total mass density in the plane 
of spiral galaxies. Van der Kruit (1981) showed that 
the width of a galaxy's gaseous disk should increase 
radially in an exponential fashion. From the thick- 
ness measurements of the gaseous disk of the edge-on 
galaxy NGC 891 by Sancisi & Allen (1979) van der 
Kruit concluded that the M / C ratio of the stellar disk 
does not change significantly with radius and that the 
dark halo can not be as flat as the stellar disk. Rupen 
(1991), using more sensitive and much higher resolu- 
tion data, found that the width of the gaseous disks 
of NGC 891 and NGC 4565 increase exponentially, 
but for NGC 891 with large scatter. Our galaxy has 
been studied in much greater detail by Knapp (1987) 
and Merrifield (1992) in HI, and by Malhotra (1994) 
in CO. These authors find that the radial scale- length 
of the midplane mass density compares well with the 
radial scale-length of the stellar disk, implying a con- 
stancy of the mass-to-light ratio for the stellar disk. 

Van der Kruit (1988) concludes that in the outer 
parts of stellar disks the local mass density is no 
longer dominated by the stellar disk so that the self- 
gravity of the gas will become important. Here I ex- 
tend van der Kruit 's pioneering work by incorporating 
the stellar and gaseous disks as well as the DM-halo 
in a self consistent way : rather than making a lo- 
cal approximation to the vertical force, is calcu- 
lated from the global mass distribution of the galaxy. 
In 



3.2 



I show that the self-gravity of the gas can 
play an important role and must be included in a self- 
consistent manner. Furthermore I will concentrate on 
the flaring behaviour in the region beyond the opti- 
cal disk where the vertical force is dominated by the 
DM-halo, and is hence sensitive to its shape. In com- 
bination with the newly achievable large sensitivities 
and the high resolution of HI synthesis observations 
(Rupen, 1991), these improved modelling methods al- 
low for the determination of the shape of dark halos 
of spiral galaxies. 

In §|2] I review some of the properties of dark ha- 



los and determine (in §2.1.) how, for a given rota- 



tion curve, the core radius and central density de- 
pend on the flattening of such an isothermal halo (see 
also Appendix A). The disk-halo conspiracy is inves- 
tigated in section §p.2.| (and Appendix B) where I 
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derive general formulas for the three unknowns (the 
disk's M / C ratio and the halo's core radius and cen- 
tral density) which determine the overall mass distri- 
bution of the galaxy. In §^ I list the assumptions 
made to determine the z-distribution of the gas from 
the galaxian potential. I compare the local and global 
approaches in ^3.1., where I also present the vertical 



force arising from the total mass distribution of the 
galaxy, Kz^tot, for several radii (an analytic solution 
to the multi-component local approximation, as pro- 
posed by Bahcall (1984), is given in Appendix C). In 
§ 3.3. I calculate how the thickness of the gas layer de- 
pends on the two free parameters {M /£ and q) and 
discuss how the thickness of the gas layer can be used 
to determine the flattening of the DM-halo. In the 
discussion, I will indicate which systems might 
be suitable for an analysis as proposed in this paper. 



2. Properties of dark halos 

From the observational fact that rotation curves of 
spiral galaxies are "flat" , it has been concluded (e.g. 
Bosma 1978; Rubin et al. 1980; Bahcall & Casertano 
1985; Begeman 1987) that there is an unseen mass 
component present in these galaxies. The standard 
assumption has been that this dark mass distribu- 
tion is spherical and isothermal, characterized by the 
core radius, i?ci ^^nd central density, phfl (equation 
I All , with q 



c/a—1). The luminous matter con- 
sists of a stellar disk, for some systems a bulge, and 
a gaseous disk. The M./ C ratios of the bulge and 
the disk are free parameters but are normally taken 
to be constant with radius. An upper limit to the 
M / C ratio of the disk, and hence a lower limit to 
the dark halo mass, can be obtained by assuming 
that the observed peak rotation velocity is due to 
the stellar disk only. A more common approach is to 
scale the mass-to-light ratios down in such a way as 
to avoid halos with hollow cores. This is commonly 
known as the "maximum-disk" hypothesis (van Al- 
bada & Sancisi, 1986). Other means of constraining 
the ratio exist : van der Kruit (1981) calcu- 

lates M-l L from the thickness of the gas layer, Efs- 
tathiou et al. (1982) invoke disk stability arguments, 
Athanassoula et al. (1987, ABP87) apply spiral insta- 
bility criteria, and Bottema (1993) uses stellar veloc- 
ity dispersion measurements. The mass-to-light ratios 
found by these authors range from 50 to 100 % of the 
maximum-disk value. 

The values for the DM-halo core radius and the 



A4 /C ratio of the stellar disk are highly correlated : 
low mass-to-light ratios require small core radii, and 
large M /£ ratios correspond to large values for the 
DM-halo core radius. In fact, many different disk-halo 
decompositions produce acceptable fits to a given ob- 
served rotation curve (van Albada et al., 1985 here- 
after ABBS85 ; Persic & Salucci, 1988 ; Lake & 
Feinswog, 1989 hereafter LF89). The dark to lumi- 
nous mass ratio may decrease with increasing mass 
(e.g. ABP87; Casertano k van Gorkom 1992; Broeils 
1992, Chapter 10). This could be inferred from galac- 
tic rotation curves : dwarf galaxies have rising rota- 
tion curves (e.g. Carignan & Freeman 1988) while the 
rotation curves of massive galaxies fall in the outer 
parts ( Casertano & van Gorkom 1991; Broeils 1992). 
For those galaxies where the rotation curves in the 
outer parts are flat, the luminous and dark matter 
"conspire" : the increase of V^^j^ and the decline in 
^disk ^^'^'^ t'^^t their sum remains approximately 
constant with galactocentric radius. 

In view of the uncertainties in the disk-halo de- 
composition, it is clear that more and a different kind 
of data are needed for a unique determination of the 
halo flattening. In section § 3.3. I show that beyond 
the optical disk the thickness of the gas layer is rather 
sensitive to the flattening of the halo. In order to in- 
vestigate this sensitivity quantitatively I construct a 
two parameter galaxy model, where these parameters 
are 7 {— the fraction of the peak rotation curve due 
to the stellar disk, sec §2.2. and Appendix B) and q, 
the flattening of the DM-halo (§2.1. and Appendix 



A). The explicit dependence of the thickness of the 
gas layer upon these two parameters is discussed in 



3.3 



2.1. Flattened dark halos 

As mentioned above, the equatorial rotation curve 
(i.e. the rotation curve in the plane of the stellar disk) 
does not constrain the actual shape of the halo. In 
Appendix A I determine a family of flattened DM- 
halo models which have the same (to within 1.4 %) 
equatorial but different polar rotation curves. Thus, 
for this family of halo models, the radial force does 
not depend significantly on the shape of the DM-halo, 
while the vertical force (and hence the thickness of 
the gas layer) does. In Figure |l|, I present the rota- 
tion curves for this DM-halo family graphically. The 
lower panel shows the rotation curves for these DM- 
halo models, while in the top panel I present the ra- 
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tio of the flattened to the round DM-halo rotation 
curve. Although the residuals show systematic be- 
havior, the amplitudes (< 1.4%) are smaller than the 
routinely obtained observational errors (e.g. Bege- 
man 1989, BE89 hereafter; Broeils 1992). In con- 
clusion : rotation curves of flattened DM-halos are 
indistinguishable from their round equivalents. This 
family of flattene d D M-hal o m odels is fully specified 
by the equations ( A8) and (A9) which relate the core 
radius and the central density to the flattening of 
the DM-halo and is graphically presented in Figure ||. 
For a given rotation curve, flattened DM-halo models 
have larger core radii and central densities than their 
round equivalents. Notice that both the central den- 
sity and the core radius have an almost linear depen- 
dence on the halo flattening for moderately flattened 
DM-halos. 

Of course the DM-halo model I have chosen to 
work with might very well be different from the true 
DM halo mass distribution. However all mass dis- 
tributions, with similar rotation curves, share the 
general feature that flatter distributions have larger 
vertical forces. Therefore, a DM halo flattening de- 
termined using the formalism outlined in this paper 
serves as an indicator of the true flattening of the 
(unknown) dark matter density distribution. 

2.2. On the disk-halo conspiracy 

As we can measure the light distribution of the 
stellar disk only, and have no a priori knowledge of the 
mass-to- light ratio of stellar disks, the relative contri- 
butions of luminous and dark matter are not known 
(ABBS85 and LF89). However, as I wiU show in Ap- 
pendix B, the galaxy mass model described above is 
fully determined by one parameter regulating the rel- 
ative importance of stars and dark matter. Following 
Bottema (1993) I choose this parameter, 7, to be the 
fraction of the peak observed rotation curve which is 
due to the stellar disk (defined by eqn. [B15]). 



Different choices of 7 result in quite different val- 
ues for the stellar mass-to-light ratio, the halo's cen- 
tral density and core radius. As an example, I present 
the dependence of these three parameters on 7 for the 
galaxy NGC 3198 (rotation curve and optical param- 
eters were taken from BE89) in Figure |^, and alge- 
braically by the equations ( pT6| ), ( pl^ ) and ( pT9| ). 
In Figure |^, I have also indicated how a 5% uncer- 
tainty in the slope of the rotation curve affects the 
results. To obtain the core radius and central density 



of a flattened dark halo, one uses the equations (B18) 
and (B19) and multiplies these values by C{q) (equa- 
tion |A8[) and H(g) (equation respectively. 

In agreement with LF89, 1 find that observed equa- 
torial rotation curves do not constrain the core radius 
of the DM-halo (or 7 in our terminology) very well. 
This is illustrated in Figure ^ where I present the 
rotation curve of a model galaxy which resembles the 
Sc galaxy, NGC 31980. Several acceptable fits, made 
with different 7's, are shown. 

In conclusion : the parameterization of the disk- 
halo conspiracy, as presented in Appendix B, provides 
a useful way to perform the disk-halo decomposition 
to acceptable accuracy, and allows for a straightfor- 
ward way to investigate the dependence of the thick- 
ness of the gas layer on the 7 value (mass-to-light 
ratio) of the stellar disk. 

3. The method 



In this section I derive the vertical distributions of 
the HI layer from the potential of the whole galaxy, 
^ galaxy {R, z) , and the equation of hydrostatic equilib- 
rium. This derivation is subject to several simplifying 
assumptions which are listed below : 

• The system is in steady state, 

• and is azimuthally symmetric. 

• The velocity dispersion tensor of the gas is sym- 
metric and round, so that 

• The equation of hydrostatic equilibrium is a 
good approximation to the vertical Jeans equa- 
tion. 

• The vertical velocity dispersion of the HI gas is 
constant (isothermal) with z- height. 

• Magnetic and cosmic ray pressures are neglected 

Obviously this is just an approximate description 
of reality. Spiral density waves violate the first two 



1 The optical parameters (L(0)=207.1 Lo/pc^ Hr = 2.3 
kpc, = 0.23 kpc), gaseous surface density distribution 
and rotation curve (¥063(2. 3/ijj) = 146.4, Vobs(10.0/iij) = 
148.0, Vgas(2.3/ifl) = 13.4 and Vga^lO.Oftji) = 38.0 
km/s) were taken from BE89 who used Kent's (1987) pho- 
tometry to calculate the rotation curve due to the stellar 
disk. The stellar disk is truncated at (van der Kruit, 
1988). 
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assumptions, while non-thermal pressure terms (mag- 
netic, cosmic ray heating, ...) are not included in the 
equation of hydrostatic equilibrium. Keeping these 
assumptions in mind, I take as the starting point the 
vertical Jeans equation in cylindrical coordinates (as 
usual, R, 6, and z denote the radial, the tangential 
and the vertical direction respectively) : 

= p{z)K.{z)~ 

to determine the vertical distribution of the gas (p(z)). 
Here Kz{z) is the gravitational force (per unit mass) 
in the -|-z-direction and cr(i?, z) the gaseous velocity 
dispersion. Due to the assumption of azimuthal sym- 
metry, the last term in equation (|^) vanishes. The 
c_Rz-term essentially measures the tilt of the velocity 
dispersion ellipsoid, for cylindrical rotation it is iden- 
tically zero. In the most extreme case the velocity 
dispersion ellipsoid points towards the galactic cen- 
ter|. 

At small zj R, ajiz is approximately {a\j^—a1z)z j R 
so that its contribution is expected to he smallQ. 
Thus, the non-diagonal terms of the velocity disper- 
sion tensor (i.e. ) will vanish. Then, the vertical 
Jeans equation reduces to the equation of hydrostatic 
equilibrium : 

J-_^ = pgas{z) K,{z). (2) 

With the assumption that the gas is isothermal in the 
z-direction it follows that : 

qalaxy 

^z.qas 1 — J ' v'^/ 

d Z d Z 

so that 

Pgas{R,z) = p,,,(i?,0)e-*— (4) 

«poexp(~zV2W^2^J. (5) 

thank the referee, Phil Maloney, for pointing this out. 
^In the Galaxy, auR and a^z, as measured by Malhotra's 
(1994) and Blitz et al. (1984) are almost equal : 7.8 ± 3 
versus 5.7 ± 1.2 km/s. 



We see that the gaseous density distribution can be 
calculated once the gaseous velocity dispersion and 
the potential of the g&\axy,<^gaiaxy{R, z), are known. 
Generally speaking, the vertical distribution of the 
gas will not be a "simple" function of z. Since the ob- 
servational data does generally not allow for a more 
sophisticated analysis than the fitting of Gaussian 
functions to the measurements, I choose to fit Gaus- 
sians to the model density distributions (eqn. j|]) as 
well^. In the remainder of this paper I will use the 
dispersion {Wgas) of this Gaussian, as defined by eqn 
1^, as a measure of the width of the gas layer. 

I incorporate three components in the present 
study : 1) a stellar disk with constant scale-height, 
2) a flattened isothermal DM-halo, and 3) a gaseous 
disk. The stellar density distribution for the stars 
must be close to the one derived from surface bright- 
ness measurements by van der Kruit & Searle ( 
vdKS81a&b and vdKS82a&b), namely : 

Ps{R,z) 

^ M/C ps(0, 0) e-^/''« sech2(z/2ze) R < Rmax 

= 0.0 R> R,na., (6) 

with Rmax = (4.7±0.7)/ifl, 2ze « /ifl,/(4.7±1.8), and 
A4 /C the average mass per unit luminosity. 



* Since the shape of the potential changes with galacto- 
centric radius, the gaseous density distribution will also 
change shape. Thus, the accuracy of the approximate 
density distribution (eqn. [^]) will also change with 
galactocentric radius. Other indicators for the width of 
the gas layer could be used as well. For example the 
normalized second moment of the density distribution, 

Wmom2 = z'^ p{z) dz I p{z) dz, is more sen- 

sitive to the high-z parts of the density distribution, yield- 
ing larger values for the "width" of the gas layer. On the 
other hand, the width could be calculated from the Full 
Width at Half Maximum {Wfwhm = FWHM/2.2,^), a 
calculation which is more sensitive to the low-z part of 
the density distribution, so that smaller values for the 
"width" are found. Significant differences between these 
three width measures indicate that the true density dis- 
tribution is not Gaussian. For the toy model discussed 
below the Wfwhm and Wmom2 width estimators can 
be as much as 20% larger and smaller than Wgas respec- 
tively (inside the optical disk). Beyond the optical disk, 
the gaseous distribution is close to Gaussian. 
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Note that Barteldrees & Dettmar 1993 (BD93) find that the optical disks are truncated at significantly smaUer 
radii : Rmax = (3.0±0.2)ft,fl, 2ze ~ /i_r/(4.0±0.3). In external galaxies it is very hard to observationally distinguish 
a sech-squared distribution from an exponential distribution (i.e. Wainscoat et al. 1989). For the Galaxy on the 
other hand, there is a large body of evidence suggesting that the vertical distribution can be better represented by 
an exponential distribution^ : Ps{z) — Ps(0) (Gilmore & Reid 1983; Pritchet 1983; van der Kruit 1986; 

Yoshi tt al. 1987). Therefore, I calculate model potentials for the exponential distribution. I chose the flattened 
isothermal distribution as proposed by SS90 (equation |A1[ |), discussed above, to represent the DM-halo mass 
distribution. 

The vertical force of a truncated stellar disk as well as of a flaring disk is calculated from the potential : 



^{R,z) = -2G / rdrp{r,0) / dz'p{r,z') / — ^ , (7) 
Jo J~oo Jo \r ~ R\ 

K,{R,z) ^ -^<^{R,z) (8) 
dz 

= -2G f7drp{r,0) Hdz' , p{r, z')iz - z')fiR, z,r, z') 

Jo i-oo ^{R + r)2 + {z- z' f ([R - r]2 + [z- z'Y) ^' 



where p{r, z') is the density distribution for which to calculate the forces. At the center of the galaxy, /(O, z, r, z') 
equals tt, for all other galactocentric radii f{R,z,r,z') is the complete elliptic integral of the second kind : / = 
E(k) = d(\)\J 1 — sin^ 0, with = 4r/[(l + r)^ + {z — z')'^]. With special care for the region around {R, z) 
this integral can be solved numerically and takes about 8 seconds per point on a SPARC 2 processor. I use an 
exponential for the vertical density distribution, p(0, z'), of the stars, while I approximate the gaseous distribution 
by a Gaussian. As neither the vertical scale- height nor the midplane density, p{r, 0), are generally analytic functions, 
I store their values in tabular form and determine the interpolated value whenever the integrating routine (QROMB 
or QROMO, Press et al. 1990) requires so. I used the double exponential stellar disk (Kuijken & Gilmore, 1989, 
hereafter KG89, their eqn. [27]) as a test case and found that the two methods of calculation agree to within 1 
part in 1,000|. 

3.0.1. The global approach 

As a first step I use the radial part of the potential to determine the structural parameters of the stellar disk 
and DM-halo by requiring that they reproduce the observed equatorial rotation curve. The radial force due to the 
gas is calculated (using the ROTMOD program in GIPSY, van der Hulst et al. 1992) assuming that it is infinitely 
thin. I use the observed photometry and a constant M./ C ratio to calculate the rotation curve due to the stars 
(also by using ROTMOD). I then perform the disk-halo decomposition outlined in Appendix B. 

I use the values for the mass-to-light ratio of the stellar disk and the core radius and central density of the 
DM-halo, as found in the disk-halo decomposition, to calculate the vertical force due to the stellar disk (equation 
j|]), and the DM-halo (SRJF94's equation [6]). These forces are integrated numerically to yield the potential, from 
which the gaseous volume density distribution, pgas{R, z), is calculated (using eqn. [Q, an assumed gaseous velocity 
dispersion, and an observed surface density distributions). 

I follow an iterative procedure, where the gaseous density distribution, calculated from an approximation to the 
true total potential, is used to estimate the gaseous contribution to the potential, which must then be added to the 

^In this paper I do not include any thick disk component. Note that these isothermal and exponential vertical distributions 
have the same slope at large z-heights. Thus for a given surface density and high-z slope, the exponential distribution has 
a twice larger value of the midplane density than the isothermal distribution (see also van der Kruit, 1988). 

^This accuracy equals ten times the accuracy obtained in the integrating routines. 
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contributions from the stellar disk and DM-halo to yield a better approximation to the true total potential, from 
which a better approximation to true gaseous density distribution can be calculated, ... etc. 
Thus, the vertical distribution of the gas is not fully consistent with the Poisson equation 

A n (R \ 1 Fr) ^ 1 dK, 

= -- + 1^9^-^^ (10) 



since I take the vertical distribution of the gas layer to be Gaussian, and because the radial force (Ffj) due to a 
thick gas layer differs from the radial force arising from an infinitely thin disk. These inconsistency are expected 
to be small because of the small scale-height and low mass of the gaseous disk. This procedure to calculate the 
thickness of the gas layer is what I term the global approach. 

In order for the global model to be self-consistent, the stellar velocity dispersion tensor (cr(i?, ^ , with ij 
all the possible coordinate combinations) has to be such that for the assumed stellar density distribution and the 
calculated total potential, the vertical Jeans equation (|^) will be satisfied. Thus, every mass model described above 
makes two predictions, firstly the radial variation of the width of the gaseous layer, and secondly the stellar velocity 
dispersion. In this paper I concentrate on the radial variation of the thickness of the gas layer. 



3.1. The Local Approach 

To date, most authors have determined the distribution of the gas above the plane from the local mass densities 
only. In this section I compare this local approach with the global approach described in the previous section. To 
separate variables I will postpone investigating the effects of the gaseous self-gravity till 



3.2 



First I describe how in the local approach the vertical force is calculated from the local mass densities and 
gradients in the rotation curve. Assuming that, in the region where the gas is found, the rotation speed does not 
depend on ^-height and that the vertical force does not depend on radius, the Poisson equation (^o|) reduces to the 
local Poisson equation: 



1 d{R Fr) dK, 



Combining the appropriate equation of hydrostatic equilibrium (pi) with the equation above (O) I find : 



^'3^' dz Jo '^"^ ■ KAttGJ \ R dR 

Setting the rotation curve term to Prot, it follows that : 



^2 



d \npgasiz) 



<^t,gas ' ' K^{z) = -AttG I dz' {pmatter{z') + Prot{z')} ■ (12) 



Negative and positive rotation densities {protj correspond to rising and falling rotation curves respectively. The 
effect of a rising rotation curve is to increase the thickness of the gaseous and stellar disks, while in a falling rotation 
curve a disk will be thinner relative to the case where prot = . Bottema et al. (1987) find that the rotation term 
can contribute significantly in the rising parts of galactic rotation curves. For the case of a flat rotation curve, 
equation (9) leads to the plane-parallel-sheet case. Equation (|l^) can be solved either iteratively or analytically (if 
more simplifications are introduced) : the multi-component and Gaussian-equivalent methods described below. 



7 



3.1.1. The Gaussian equivalent method 

In order to arrive at the (analytic) Gaussian- 
equivalent method, we simplify matters further by 
assuming that the self-gravity of the gas is negligible 
and that the total mass density does not vary signif- 
icantly with height above the plane. I then arrive at 



d hipgas{z) 



z,gas 



-Att G ptotiO) z, 



so that the vertical distribution is Gaussian with dis- 
persion Wgas,G ■ 

Pgas{z) = Pgas(0)exp(-zV2W^2 .g). 



For which the following relation holds 



gas,G 



' z,gas 



TT G ptoM 



(13) 



In the Gaussian-equivalent method one then assumes 
that the Gaussian width {Wgas^c) can be equated 
with the observed width Wgasi^ FWHM/2.35, as 
defined by eqn. [^), so that the inferred midplane 
mass density (pm/) is given by : 



Pinf 



AttG 



^ z,gas 
* ^ gas 



(14) 



Inside the stellar disk, the Gaussian width (Wgas.c) 
will always be smaller then the true width because, for 
any density distribution, the volume density decreases 
with z-height. The inferred midplane mass density, as 
calculated from the true width Wgas (which is larger 
than Wgas,G) and equation (Q, will therefore under- 
estimate the true midplane mass density. This effects 
is illustrated in Figure ^, where I plot the ratio of the 
inferred midplane mass density (pm/) and the actual 
midplane density (ptot)- For several mass models re- 
sembling NGC 3198, this method recovers 50 to 110 
% of the true midplane density. Especially for the 
round halo case {q = 1.0) we clearly see the effects of 
the truncation of the stellar disk at i? = 6/i/j. The 
truncation of the stellar disk of NGC 3198 occurs at 
a relatively large distance {6hn versus on av- 

erage, vdKS81a&b and vdKS82a&b, and BD93) so 
that it's effect may be relatively small for NGC 3198. 
On the other hand, the effects of the truncation of 
the stellar disk are smaller if the vertical distribu- 
tion of the stars is sech^ rather than exponential (i.e. 
less concentrated towards the midplane). It should 
be noted that these results depend somewhat on the 



width estimator being used (see footnote for ex- 
ample the radial variation of pi„f / Ptot arrises partly 
because the vertical distribution of the gas gradually 
changes from a distribution which is strongly centrally 
peaked (in the inner two scale-lengths) to a Gaussian 
(beyond the optical disk). 

Including the rotation density will improve the 
Gaussian-equivalent method beyond the truncation of 
the optical disk (right-hand side panel). The radial 
variation in the density ratio beyond the optical disk 
results from small scale structure in the model rota- 
tion curve. However, this is possible only if the sum 
of the rotation and DM densities is larger than zero 
(see eqn. and Appendix ^]). For example, this is 
not the case for a maximum-disk model of NGC 3198 
in the inner two radial scale- lengths. 

Although its accuracy varies systematically with 
radius, the Gaussian-equivalent method works reason- 
ably well, especially beyond the stellar disk. 

3.1.2. The plane parallel sheet case 

The difference between the plane-parallel-sheet case 
and the global approach is most easily illustrated by 
plotting the vertical forces for the two cases. In Fig- 
ure H I plot the Kz{z) forces, as calculated for a 
NGC 3198-model, for a selection of radii. The dashed 
line is the vertical force due to a plane parallel sheet 
of stars with an exponential vertical distribution and 
a column density appropriate for the indicated radius. 
The open squares represent the force generated by a 
radially truncated exponential disk. The full line rep- 
resents the Kz (z) generated by the combined density 
distribution of the stellar disk and the DM-halo. For 
the combination of parameters, maximum-disk and 
round DM-halo (top-left panel), the plane-parallel- 
sheet case overestimates the vertical force in the in- 
ner two scale-lengths of the stellar disk, because the 
plane parallel sheet is a bad approximation to the 
fast radially declining surface density (there is very 
little surface area at the central surface density) . Fur- 
thermore, the large negative prot associated with the 
steeply rising rotation curve (Figure ^ decreases the 
effective volume density. Beyond about 2 hn the po- 
tential due to the interior stellar disk becomes more 
important and will eventually dominate the potential 
generated by the local column of stars. Additional 
model calculations (top-right panel : 7 = 0.9, q — 0.1 
; bottom-left panel : 7 — 0.6, q = 1.0 ; bottom-right 
panel : 7 = 0.6, q ~ 0.1) show that the stellar disk has 
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a dominating contribution to the vertical force (over 
the extent of the stellar disk) if and only if the stel- 
lar disk is close to "maximum" and if the DM-halo is 
close to round (see also §B.4 and Figure p^ . 

3.1.3. Inferences 

From the figures |[ || and |^ we learn the following 



— At large galactocentric radii (i? > Ghu), is al- 
most linear for the range of models considered here, 
resulting in an almost Gaussian vertical distribution 
for the gas. Notwithstanding the low stellar surface 
densities at these large distances, the contribution of 
the stellar disk to the total vertical force can be sub- 
stantial : even at eleven radial scale-lengths the stellar 
disk contributes about 20 % of the total vertical force 
(for the round DM-halo cases). The self-gravity of the 
gas (s ee Figure 0) , which is discussed in more detail in 
§3.2., can be important (up to 50% of the total force) 
in this region as well. 

— In the transition region (4/ijj < R < Ghjf ) the con- 
tribution of the stellar disk, the gaseous disks, and the 
dark halo to the total vertical force can be of com- 
parable magnitude (for the mass model with maxi- 
mum disk, round halo, and typical gaseous surface 
densities. See Figure Stellar disks are truncated 
(vdKS81a&b; vdKS82a&b; BD93) causing a sharp 
discontinuity in the stellar densities. However the po- 
tential as generated by the total stellar mass distri- 
bution will change much less abruptly. Thus the local 
approaches will strongly overpredict the thickness of 
the gas layer in this region (see Figure Moving 
away from this corner (7 = 0.9, q = 1.0) in parame- 
ter space increases the halo's contribution to the total 
potential significantly. 

— Inside the optical disk, the potential is dominated 
by the stellar disk only in case the disk is close 
to "maximum" and the halo is almost round. For 
the other extreme case (minimum-disk-|-flat-halo) , the 
potential is completely dominated by the DM-halo. In 
the two remaining corners of parameter space, stellar 
disk and DM-halo contribute about equally to K^- In 
all of these cases, the density distribution of a tracer 
population will be almost Gaussian as long as the ve- 
locity dispersion of the tracer is so low that most of 
the tracer resides at z-heights below the "knee" in Kz 
(at two to three stellar scale- heights) . For a tracer 
population like the atomic gas with a velocity disper- 
sion of roughly 7 km/s, the entire extent of the optical 
disk falls in this regime (see Figure Kw. Isothermal 



tracers with a velocity dispersion larger than two or 
three times (T^^gas can move beyond the knee in the Kz 
curve so that their density distribution can no longer 
be adequately described by a single "simple" function 
(cf. the thick disk of the Milky Way [e.g. Gilmore & 
Reid 1983]). 

— As expected, both the amplitude and the shape 
of the vertical force are almost independent of the 
7 value of the stellar disk if the DM-halo is highly 
flattened. In case of a round DM-halo, different 7 
values for the stellar disk do result in different shapes 
for the Kz curves. 

As the amplitude of the knee in Kz depends lin- 
early on central surface density (quadratically on 
Vrot.max, equ. |B10|), the thickness of the gas layer 
will be smaller in high rotation curve galaxies and 
larger in dwarfs (see also Appendix D. ). 

Two radial regimes are poorly modeled by the 
Gaussian-equivalent as well as the plane-parallel-sheet 
local approaches : the region where the rotation curve 
rises so sharply that {prot + Pdm < 0), and the region 
where the stellar disk truncates. Special treatment 
of these regions can be avoided by using the global 
approach. For those regions where the rotation curve 
is almost flat, which are far away from the truncation 
of the stellar disk, the Gaussian-equivalent method 
works reasonably well. 

3.2. The self-gravity of the gas 

Rupcn (1991) flnds in NGC 891 as weU as in NGC 
4565, that the HI volume densities become compa- 
rable to the stellar volume densities beyond about 2 
radial scale-lengths. For the Milky Way this is true 
at the solar circle. The self-consistent solution for the 



multi-component method as found in Appendix B.4, 
(eqn. [ C28 ) indicates that the importance of the self- 



gravity of the gas scales with the ratio of gaseous to 
stellar kinetic energy. At the solar circle this ratio 
is about 0.1-0.2, so that the self-gravity of the gas is 
expected to become increasingly important at larger 
radii. 

The importance of the self-gravity of the gas is 
illustrated in Figure 0, where I plot the K^ due to 
the total (stars -|-DM-|-gas) potential (the full line), 
and the forces arising from the global DM, stellar and 
gaseous density distributions (dotted, short dashed, 
and long dashed lines respectively). The self-gravity 
of the gas has be en calculated using the procedure 
outlined in §3.0.1. 
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3.2.1. The multi- component method 

In case prot and the self-gravity of the gas are 
neglected it is rather straightforward to show (van 
der Kruit, 1988) that different isothermal components 
subject to the same cylindrically symmetric gravita- 
tional field have vertical distributions which are pow- 
ers of each other. Including p^ot, the self-gravity of 
the gas, and the dark matter halo leads to the multi- 
component method, first used by Bahcall (1984) to 
determine the amount of disk-like dark matter in the 
solar neighborhood. In Appendix [c], I present an an- 
alytic solution (eqn. | C28 |, which needs to be solved 
iteratively) to equation (p_2[) . In the same Appendix I 
also show that in the multi-component method the 
distribution of the gas is a power of the distribu- 
tion of the stars as well. In Figure ^ I compare 
the model widths calculated in the multi-component 
method (where the vertical distribution of the stel- 
lar disk is fixed to a sech-squared , see Appendix 
PI with those calculated using the global approach. 
The top plot shows the percentage difference between 
the global and local approaches. The structure in 
these difference curves mostly arrises due to small 
scale structure in the observed rotation curve. Thus, 
application of the multi-component method could be 
misleading. For example, if the flaring measurements 
extend to a radius at which the local approach pre- 
dicts a local maximum in the observable width, too 
round a DM-halo would be inferred. 

I used the multi-component method to compare 
models of NGC 3918 with a varying amount of gaseous 
self-gravity. The results are presented in Figure H. In 
the left-hand set of panels I show the width, Wgas (~ 
FWHM/2.'ib, as defined by eqn. [§), of the gas layer 
as a function of galactocentric radius for three differ- 
ent halo flattenings (bottom panel, q = 1.0; middle 
panel, q = 0.5; upper panel, q — 0.1). Four different 
curves are presented, the different lines corresponds 
to differently scaled gaseous surface densities : 1/100, 
1, 5, and 10 times the measured surface density for 
the full line, the dotted line, the short dashed line and 
the long dashed line respectively. The models are cal- 
culated!^ for a gaseous velocity dispersion of 7.5 km/s 
and a stellar disk with a 7 value of 0.93 (i.e. about 
maximum-disk). For the scaled gaseous disks I plot 
(right-hand set of panels) the ratio of thickness of the 



^The effect of the scaling of the gaseous surface density 
has been taken into account while performing the disk- 
halo decomposition. 



scaled disk and the thickness of the no-self-gravity- 
disk. From this plot we see that not including the 
gaseous self-gravity will result in an overestimate of 
the width of the gas layer (by approximately 10% to 
20%). For the predicted gas-layer width to equal the 
observed widths a more massive DM-halo is required 
(eqn. |S9|), so that neglecting the self-gravity of the 
gas will result (using eqn. p^) in a smaller value (by 
approximately 20% to 40%) for the halo flattening, q. 
Thus, in the local approach as well as in the global 
approach (see Figure |l^) , the self-gravity of the gas 
is found to have a comparable effect on the inferred 
DM-halo flattening. 

From these plots I also conclude that for galax- 
ies with relatively large gaseous surface densities the 
thickness of the gas-layer is not sensitive to the flat- 
tening of the dark matter at all. Thus, the self-gravity 
of the gas may be an unimportant or a significant con- 
tributor to the vertical force (for early and late type 
systems respectively). 

3.3. Halo shape and gas layer thickness 

In order to investigate the dependence of thickness 
of the gas layer upon the flattening of the DM-halo 
and the M-j L ratio of the stellar disk, I determined 
Pgas{R,z) (and Wgas{R) thcrcfrom) for a range of 
galaxy models. 

In Figure |l^ I plot the width, 
Wgasi- FWHM/2.35, as defined by eqn. [|), of 
the gas layer in a 7 = 0.9 (i.e. ~ maximum-disk 
) disk-halo combination for several halo flattenings. 
The left-hand panel was calculated for a double expo- 
nential disk. In the middle panel I incorporated the 
truncation of the stellar disk (at R/hu = 6), while 
the self-gravity of the gas is included in the right- 
hand panel. The effect of the self-gravity of the gas 
is different for the different disk-halo models, and is 
strongest between 4 and 8 radial scale-lengths : be- 
tween 20% and 40% for the NGC 3198 models. Be- 
yond several radial scale-lengths the width of the gas 
layer depends strongly (roughly proportional to y^) 
on the flattening of the DM-halo. Inside the stellar 
disk, the width of the gas layer is largely independent 
of the DM-halo flattening because the potential there 
is dominated by the stellar disk. However, for mod- 
els with flatter halos and/or lower-7 stellar disks, the 
halo can be of importance inside the stellar disk (see 
§3.2.1., and Figure It is important to keep in mind 
that the self-gravity of the gas is more important for 
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the rounder DM-halo models (§ p.2j and compare the 
Figures |^a and ^b) . So the "signature" of the halo 
flattening is somewhat reduced if the self-gravity of 
the gas is important (see also Figure ||). 

Is the thickness of the gaseous layer dependent on 
the mass-to-light ratio of the disk? In order to answer 
that question, I performed the calculations described 
above for several 7's. The results (for q = 1.0) are 
presented in Figure [l[ We see that, as expected, 
inside the stellar disk a gas-layer immersed in a more 
massive stellar disk will be comparatively thinner. At 
large galactocentric radii the thickness of the gas layer 
becomes virtually independent of the M./ C ratio of 
the stellar disk. The same conclusion can be reached 
from Figure]^, where we see that, beyond the optical 
disk, the derived midplane mass densities for different 
M./C -disks are nearly equal. 

3.4. Effects of incomplete modelling 

It is instructive to see how uncertainties in gaseous 
velocity dispersion and uncertainties in the derivation 
of the vertical force affect our estimates of the DM- 
halo flattening. Another point of interest is to inves- 
tigate what differences may arise when using the local 
instead of the global approach. 

To summarize the following sections : the errors 
introduced by using the Gaussian-equivalent method 
(typically a 10-30% too round halo) tend to cancel 
out those introduced by neglecting the self-gravity of 
the gas (20-60% too flat a halo). 

3.4-1- Uncertainties in the velocity dispersion 

In many practical situations it may be that mea- 
surements of the gaseous velocity dispersion are not 
available. This severely limits our ability to deter- 
mine the flattening of the DM-halo. In some face- 
on galaxies the azimuthally averaged gaseous velocity 
dispersion is observed to decline radially from around 
12 km/s in the center to remain constant at 7 km/s 
beyond i?25 (van der Kruit & Shostak 1982; Dickey 
et al. 1990). In other systems, the region beyond i?25, 
as well as the inter-arm regions appear to have a con- 
stant velocity dispersion of 7 km/s, while on the arms 
values around 13 km/s are reported, ( Shostak & van 
der Kruit 1984; see Kamphuis 1993 for a review). In 
summary, the measurements to date indicate that the 
gaseous velocity dispersion beyond the optical disk 
has a value of around 7 km/s. 



Overestimating the gaseous velocity dispersion by 
a factor of (1 -|- x) causes an overestimate of the mid- 
plane mass density by a factor of (H-x)^. Beyond the 
optical disk, too flat a DM-halo distribution would be 
inferred (by the same factor (1 -f x)^ ). 

3-4-2- Uncertainties in the vertical force 

Here I estimate the effects of neglecting certain 
components to the vertical force (e.g. the self-gravity 
of the atomic gas, magnetic forces, or the gravita- 
tional force due to some remote point source). A 
mass-model which uses an underestimate of the true 
vertical force requires a larger DM density and thus 
a smaller values of the halo flattening to reproduce 
the observed gas layer widths. Thus, underestimat- 
ing the vertical force by a factor (1 -t- y) results in too 
small an estimate for the halo flattening (by a factor 
[1 + y]). The inverse is true when working with an 
overestimate of the vertical force. 

Neglecting the self-gravity of the neutral atomic 
gas, which c an t ypically contribute 30 % to the ver- 
tical force (§ 3.2. and Figure 0), will result in a 30 % 
too small inferred halo flattening. Since the known 
molecular gas content of galaxies rarely exceeds 10% 
of the stellar mass (which contributes roughly 20% to 
Kz at large distances, see §3T. and Figure ^), and 
is confined to the inner parts of the galaxy like the 
stars, I expect that neglecting the molecular gas will 
result in an underestimate of q by at most 2%. 

3-4-3- Using the Gaussian- equivalent method 

Application of the Gaussian-equivalent method to 
an observed set of gas layer widths underestimates 
the inferred midplane density by a radially varying 
amount (see §3T^, and Figure As a result the 
inferred halos will be too round and have core radii 
that are too large. 

4. Discussion and Conclusions 

In the previous sections we have seen that a flat- 
tened dark halo imprints a very specific signature on 
the flaring of the gas-layer (Figure [l0|). Measurement 
of this flaring allows for the determination of the DM- 
halo flattening. 

Although the global approach is limited by some of 
the same assumptions as the local approach (isother- 
mality of the gas, zero non-thermal pressure terms, 
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constant M-j L -ratio of the stellar disk, alignment be- 
tween DM-halo and stellar disk, ...), it treats some 
other aspects much better (radial gradients in stel- 
lar surface density, the effects of the rising rotation 
curve, the non-cylindrical symmetry of the potential, 
...). In particular the region where the rotation curve 
rises sharply and the region where the stellar disk 
truncates are better analyzed using the global ap- 
proach. In those parts where the rotation curve is 
'flat' and which are not too close to the truncation 
of the stellar disk the Gaussian-equivalent and multi- 
component methods can be used to reasonable accu- 
racy. However, the global approach is less sensitive 
to small scale irregularities (Figure To obtain re- 
liable results, any applied method must include the 
self-gravity of the gas. 

Therefor e, I p ropose that the more robust global 
approach (§3.0.1.) should be used in programs aimed 



at determining the flattening of the dark matter halo. 
For a set of global mass distributions with varying 
halo flattening one can calculate the flaring behaviour 
of the gaseous disk. The DM-halo model exhibiting 
the best agreement between calculated and observed 
gas-layer widths then yields the halo flattening. The 
accuracy of the determined flattening can be gauged 
by varying the structural parameters of the individual 
mass components. Since the velocity dispersion of the 
gas is a crucial parameter in the analysis, care must 
be taken to determine its radial dependence. This 
will be done done in a forthcoming paper where I will 
analyze the flaring gas layer of thee almost edge-on 
galaxy NGC 4244. 

4.1. Suitable systems 

Nearby edge-on spiral galaxies with large (extend- 
ing beyond the optical disk) neutral hydrogen en- 
velopes are possible targets for studies of this kind. 
Symmetric systems are preferred since it is more likely 
that such systems are in equilibrium. Furthermore 
suitable systems should not exhibit strong warping. 
For systems with moderate warps, the kinematical 
information can be used to disentangle flaring from 
warping. Unfortunately these kind of systems might 
be rare. 

In order to investigate the possibility of determin- 
ing the thickness of the gas-layer in less inclined sys- 
tems, I simulated HI observations of two models sim- 
ilar to NGC 3198 (inclined by 72° w.r.t. the hnc of 
sight). In the thin disk model, I constrained the gas- 



layer to have a FWHM of 100 pc everywhere, while 
in the flaring disk model the gas-layer has a verti- 
cal distribution calculated using the multi-component 
method (including the self-gravity of the gas, eqn. 
]C28| ). Both models were calculated with a gaseous 
velocity dispersion of 7.5 km/s. For a given set of 
observational parameters (15" beam 700 pc] and 
2.56 km/s velocity resolution) I calculate the intensi- 
ties originating from all points in the model galaxy for 
all velocities. Taking the normalized second moment 
with respect to radial velocity yields the apparent ve- 
locity dispersion for all points in the model galaxy. 
In Figure |l^ I present the two model velocity disper- 
sion maps. The systematic differences, of the order 
of 5 km/s, arise due to the fact that a line of sight 
through a flaring disk samples a wider range in posi- 
tion angle and galactocentric radius than the line of 
sight through a thin disk. In principle, the gaseous 
velocity dispersion can be determined after the layer- 
thickness-efFects are corrected for. The thickness of 
the gas layer can be derived by comparing individual 
channel maps^. Contour maps of a particular chan- 



nel are presented in Figure 13 : in the top panel the 



thin disk map is shown, the middle panel shows the 
thick disk equivalent while the difference between the 
two is contoured in the lower panel. Due to its larger 
vertical scale-height, the projected width of the thick- 
disk channel map should be larger than that of the 
thin-disk, while the average surface brightness will be 
reduced. This signature of the thick disk is easily rec- 
ognized in the displayed contour maps. Furthermore, 
contrary to edge-on galaxies, the flaring can be mea- 
sured in many different channel maps so that local 
irregularities (e.g. spiral streaming motions, bubbles, 
...) are less likely to hamper the determination of 
the radial behaviour of the flaring (cf. Oiling & van 
Gorkom, 1995). 

Therefore I expect that it must be possible to ex- 
perimentally determine both the gaseous velocity dis- 
persion and the thickness of the gas-layer for moder- 
ately inclined systems as well. 

4.2. A control sample 

It might be possible to test the importance of ne- 
glected physics in galaxies where the rotation curve 
falls in Keplerian manner. In such systems the con- 



*A channel map is an image (of a galaxy) taken in a "nar- 
row band filter" : 2.56 km/s wide in this case. 
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tribution of the dark matter to the global potential is 
presumably small, thereby eliminating the gas layer 
thickness dependence on the shape of the DM-halo. 
As a result the thickness of the gas layer should de- 
pend on the total visible mass of the galaxy and 
the gaseous velocity dispersion only. Some galaxies 
exhibit falling rotation curves, albeit not in Keple- 
rian manner (Casertano & van Gorkom 1991; Broeils 
1992). 

The measured width of the gas layer can then be 
compared with model predictions. If the gas layer 
is observed to be thicker than predicted, one might 
conclude that magnetic and cosmic ray pressures are 
important. On the other hand, a gas layer thin- 
ner than predicted indicates that either the ionizing 
extra-galactic radiation field is important (e.g. Mal- 
oney 1993), or that an unseen disk-like DM compo- 
nent is present (i.e. the "clumpuscules" proposed 
by Pfenniger et al. 1994a, and Pfenniger & Combes 
1994b). 
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4.3. Summary 

I have investigated the use of accurate determi- 
nation of the thickness of the Hl-layer in external 
galaxies. In combination with constraints set by the 
gaseous velocity dispersion and the rotation curve, 
I find that in principle it is possible to determine 
the shape of DM-halos. The gaseous self-gravity con- 
tributes significantly to the vertical force field (up to 
40%) but can be modeled accurately using the global 



approach (§3.0.1.) outlined in this paper. The thus 
determined halo flattening is independent of the par- 
ticular choice of disk-halo decomposition. 

In practice (not necessarily edge-on) systems with 
extended Hl-disks and regular velocity fields should 
be used, provided that they are observed with high 
angular and velocity resolution. 
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Appendices 

A. Rotation curves of flattened DM-halos 

As mentioned before, it is not possible to constrain the DM-halo flattening by measurements of the 
equatorial rotation curve since both round and flattened mass distributions can generate the same rotation 
curve. In this Appendix I show that the DM-halo model used by Sackett &: Sparke (1990) : 



(where Rc, Ph,o and q (= c/a) are the DM-halo core radius, central density and flattening respectively) 
defines a one-parameter family of halo models with essentially the same equatorial rotation curve if core 
radius and central density depend in a specific way on this one parameter, q (equations [A8| and [ AE] 
respectively). 

SRJF94 give a general relation for both the vertical and the radial force due to such a dark halo (their 
equation [6]). Evaluating this relation in the midplane, I find that the equatorial rotation curve arising 
from a flattened DM-halo is given by : 

yhaio{R'^Phfl,Rc,q) = Vi^aio{Qo;ph,o,Rc,q) X Q{R;Rc,q) , (A2) 

with 



Q{R] Rc, 



( M \( ^„ot.nf^'+''-/'"-W)|. (A3) 



[ \&ict&ng[q)J \R'^ + q^)Rl{q) ) \ q^Rl[q) 

and 

VLio{^;phfi,Rc.q) = 4irG ph,o{q) Rciq) f{q) , (A4) 
while f{q) and g{q) are given by : 



M arccosg\ / v^l - 

= [-7T=W) = [-^) ■ 

It is easily verified that equation ( [A2[) indeed yields the familiar 

Vj^{R) = 14^(00) {1 — (Rc/R) arctan (R/Rc)} , the rotation curve due to a round halo. 

Changing the flattening of the halo should not change the halo rotation curve significantly, that is to 
say : Vi^^i^{R;q) VhaioiR'^ 1) radii. Noting that /(I) = 1.0, I find that at infinity the following 

relation holds : 



(A6) 
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Requiring that the round-halo rotation curve equals the flattened-halo rotation curve at R = Rc{q) leads 
to the following relation : 





Noting that the RHS of this equation varies between 1.28 and 1.00 for q e [0.1, 1.0], so that Rc{q)/Rc{l) 
is constrained to the interval [1.00, 1.242], I expand (to second order) the left-hand and right-hand sides of 
this equation with respect to Rc{q)/Rc{l) and q respectively. The resulting quadratic equation is solved 
to obtain the core radius of the DM-halo as a function of flattening. The central density is then found 
from equation ( |A6D . After substituting the solutions for Rc{q) and Phfi{q) in the equation for the rotation 
curve (equation |A2| ) it is found that the rotation curves of flattened DM-halos reproduce the rotation 
curve of a round halo to within 7 % (for q =0.2). This error is maximal for i? = 0, changes sign at 
R = Rc{q) and slowly drops from -2% at two core radii to 0% at infinity. Since these deviations are large 
in the region where rotation curves are determined, such accuracy is not acceptable for our purpose. A 
little experimentation leads to a solution which is more acceptable : our best approximation ( errors < 
1.4%) to the dependence of halo core radius and central density on its flattening q is graphically presented 
in Figure and algebraically below : 



and 



R^(q) = i?c(l) X (l.209 - 0.332g + 0.123g2 

= Rcil) X C((?) , (A8) 

Ph,o{Q) = Ph,o{^) X ( /(g)c(g)2 ) + O-O^^^'^ - 0-0073g2) 

= P.,o(l) X nq) ^ ( , (A9) 



where limg^i Rc{q) = Rc{^) and lim^^i Ph,o{q) = Ph,o{^)- Substituting ([A9[) , (|A8D and (|A5|) into equation 



([Al|), I find that at large distances the local DM-halo densities are roughly proportional to 1/ f{q) ~ 1/q. 

B. On the Disk-Halo conspiracy 

In order to calculate the potential of the galaxy, the mass distributions of both the stellar disk and 
the DM-halo need to be parameterized in a convenient way. In this section I describe how the DM-halo 
structural parameters are related to those of the stellar disk. Given the photometry of the stellar disk 
and the observed rotation curve, a disk-halo mass model is fully determined by one parameter which 
regulates the relative importance of stars and dark matter. Following Bottema (1993) I choose the 
parameter, 7, to be the fraction of the peak observed rotation curve which is due to the stellar disk 
(defined by eqn. [|B15|]). With such a parameterization of the disk-halo decomposition available, it is 
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possible to investigate the effects of the stellar mass-to-light ratio upon the flaring of the gas layer. In 
accordance with van ABBS85 and LF89, I find that a wide range of stellar mass-to-light ratios (7's) 
produce acceptable "fits" to observed rotation curves. 

Since the light distribution of the stellar disk is relatively well known, the first step in tackling the 
disk-halo conspiracy is the description of the rotation curve due to the stellar disk. When assuming a 
constant mass-to-light ratio, the stellar density distribution can be described, to first order, as radially 
exponential and vertically thin. The amplitude of the rotation curve due to such a thin exponential stellar 
disk depends only on its central surface density and radial scale length. The shape is fully determined by 
the radial density distribution. Casertano (1983) has shown that both the amplitude and shape of stellar 
rotation curves are sensitive to the thickness of the stellar disk as well as to the truncation radius, Rmax 
(vdKS81a&b; vdKS82a&b; BD93). Therefore, the disk-halo conspiracy too will vary with the details 
of the stellar mass distribution. While the amplitude of the inner stellar rotation curve depends rather 
strongly on the thickness of the stellar disk, the outer rotation curve is virtually independent thereof. 
I will now investigate the dependence of the stellar rotation curve on the structural parameters of the 
stellar disk {Hr, Rmax, and Ze). 

I calculated stellar rotation curves, using the program ROTMOD in GIPSY (van dcr Hulst et al. 1992), 
for several truncation radii and disk thicknesses, and vertical distributions. As it turned out, the rotation 
curves of such stellar disks can be parameterized by a rotation speed at an inner (Rinn) £^nd an outer 
{Rout = X Rinn) radius : 



^rot.stars {Rinn) = fm7rGhRL{0)M//: (BIO) 
Vrot,starsi,Rout) ~ fx Vrot,stars (Rinn) , (Bll) 

where L(0) is the central surface brightness, Ai/C the mass-to-light ratio, and hji the radial scale length of 
the stellar disk. I take Rinn = 2.3/tR, because a typical stellar disk (truncated at 4.5/tR and hn/ze = 0.1) 
has its maximum rotation speed at this distance. The constants and fx depend upon the the vertical 
distribution and the radial truncation of the stellar disk. For the outer radius I used three different values, 
Q.OhR, 8.0/iR and W.OhR. The results are presented in Table 1 below for a range in parameters covering 
the observed range in thickness and truncation radius ( vdKSSla&b; vdKS82a&b; Bottema et al. 1987; 
BD93). 

The listed values are the average of the rotation curves for a vertically exponential and a vertically sech^ 
stellar disk. With this choice for the two vertical distributions, sech-squared and exponential, I bracket 
the suggested range of suggested distributions (van der Kruit, 1988). The uncertainty concerning the 
true vertical distribution translates directly into uncertainties, A/, in and fx- Experimentally I find, 
Afm « (2A/6) « (2A/8) w (2A/10) = 0.6%, 1.2%, 1.7% and 2.3% for zJhR = 0.05, 0.10, 0.15 and 0.20 
respectively), so that /^^^ = f + Af and f^^*^^"^ = f — Af. Thus, if the structural parameters of the stellar 
disk are well known the errors on the /-values are small. In situations where the thickness of the stellar 
disk are not known (as for all non edge-on systems) one must use the average values, fm and fx, which 
arc listed in the columns labeled ^^Zg/hR = averaged^ The uncertainties on these averages are typically 
4.5, 2.0, 2.0 and 2.0 %. If also no information is available on the location of the truncation radius, then 
fm, fe, fs, and /lo should be representative of the range in allowed values : I suggest fm = 0.73 ± 6%, 
fa = 0.68 ± 10%, fs = 0.57 ± 10% and /lo = 0.50 ± 10%. 
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B.l. The mass-to-light ratio of the stellar disk 



Now I turn to the dark halo contribution to the observed rotation curve. As usual, the (square of the) 
observed rotation curve is generated by the quadratic sum of the constituting mass components : 



^rot,obs{^) — ^rot,stars{^) + K^ot,buigei-^) + ^rot,gas{^) + Kot,ha/o(-^) • (^12) 



Now I define the gas-bulge-corrected observed rotation speed as the square root of (|B13[) . 



so that : 

^rot,obsi^) = ^rot,stars{^) + ^rotMlo^^) ■ (1^14) 

Furthermore I define 7 and /3x(= the slope of the rotation curve) : 

^ _ Vrot,stars{2.3hli) ^ _ Vobsjx hjj) (B15) 

Vrot,obs(2.3/lij) Vofes(2.3/lii;) 

Given the observed radial scale length and central surface brightness of the stellar disk, the mass-to-light- 
ratio is a function of the stellar contribution to the peak rotation curve (7) only and is given by (from 
the equations pToj and pT5[] ) : 



\Jm TtG hR L[U) J \ Jm,0.722/ifl,l-^100(Uj J 

where Vobs,ioo(2-3/ii?) is the gas-bulge-corrected observed rotation speed at 2.3 radial scale-lengths in 
units of 100 km/s, Lioo(O) is the central surface brightness of the stellar disk in units of 100 Lo/pc^, /i_r,i 
is the radial scale-length in units of 1 kpc and /m, 0.722 in units of 0.722 (the /m-value appropriate for a 
typical stellar disk with hn/ze = 10.0 and Rmax/hn = 4.5). 

B.2. The core radius of the DM-halo 

For a given observed rotation curve and a calculated stellar rotation curve shape, the core radius is 
determined uniquely by 7. Applying equation ( pi4| ) at 2.3 and x Hr, eliminating the terms depending 
on the stellar rotation curve, and using equation (|A2D with q = 1,1 find that the following relations hold 



iP'. - if. if 
1-72 



(B17) 



17 



The approximation ( [B17| ) is valid for small ranges in Rc/hn only. I list these ranges, together with 
the ranges for the left-hand side of ( B17| ) and the appropriate a^-values in Table 2. These piecewise 
approximations recover the right-hand side of equation (B17) to within 0.2 %. No effort has been made 
to match the pieces smoothly. 

For any assumed 7 and measured (5x and fx one can calculate the left-hand side of equation (|B17 ) and 
find the appropriate a^-values in Table 2 from which TZc = Rc/hR follows : 



t;^){t,P.J.) « ^ ^ ^ ^ (B18) 

B.3. The central density of the DM-halo 

The DM-halo central density is calculated from ( |B14[ ), ( pl5| ) and ( pT8| ). With 



it follows that 



Ph,o{r,Px,L) = (0.185 Mo/pc 



^3 X / Vobs, 100 {2. 3hji) 



(1-72)1 fT^ch)\_,_f 2.3 



7^c(7)2 II V 2.3 ; V7^c(7) 



arctan — — \ , (B19) 



where TZc, the normalized core radius of the DM-halo, depends strongly upon the choice for 7, and only 
weekly upon the shape of the rotation curve (via (3x) and the shape of the stellar mass distribution (via 

As an example I present the dependence of the mass-to-light ratio of the stellar disk (lower right panel 
of Figure ^), the central density (upper left panel) and core radius (lower right panel) of the DM-halo 
on 7 for NGC 3198 (photometry and rotation curve taken from BE89). The error-bars were calculated 
assuming the most common observational status : the truncation radius of the stellar disk is known but 
the thickness of the stellar disk is not. In this case, the accuracy with which one can determine the 
mass-to-light ratio of the stellar disk and the central density and core radius of the DM-halo are typically 
3, 0.5 and 4 % respectively. If the thickness of the stellar disk is also known (as would be the case for 
edge-on galaxies), the errors decrease substantially, while they blow up if neither the truncation radius 
nor thickness of the stellar disk are known. For galaxies which are about "maximum-disk", and have 
more or less flat rotation curves. Figure |^ shows that the core radius of the DM-halo will be 7 it (a factor 
of two) times the radial scale-length of the stellar disk, not too different from the maximum-disk-fits 
performed on real galaxies (Broeils 1992). Note that while the core radius of the dark halo approaches 
infinity as 7 nears unity , the dark to luminous mass (upper right panel) is not singular. 
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I tested the procedure outlined above on simulated data, where I added (in quadrature) a stellar and 
a DM-halo rotation curve to represent a simulated total rotation curve. For falling, flat, as well as for 
rising simulated total rotation curves I was able to accurately recover the input parameters, provided 
that the slope of the total rotation curve is determined over as large a range in radius as possible. 

As already pointed out by ABBS85 and LF89, mass models with a range (a factor of ten in some cases) 
in stellar mass-to-light ratios are consistent with observed rotation curves. I find, attempting to fit model 
rotation curves with "faulty" 7 values (or Ai /C -ratios) that a wide range of stellar mass-to-light ratios 
is possible if small (~ 2 %) differences between input and recovered rotation curves are allowed to occur. 
The differences between input and recovered rotation curves increase significantly inside the inner point 
and beyond the outer point which were used for to determine the slope of the rotation curve. This is 
illustrated in Figure ^ where I show the resultant fits to the observed rotation curve of NGC3198 (BE89) 
for several 7's. I took the outer radius at 10 radial scale lengths. Although all 7 values between 0.93 
and 0.53 produce similarly acceptable fits, the low-7-disks proposed by Bottema (1993) "fit" somewhat 
better. Since the velocity field of the inner parts of NGC 3198 deviates from circular rotation (Corradi 
et al. , 1991), the discrepancies between observed and model curves (in the inner parts) can not be used 
as an indication of goodness of fit. 

Thus the method presented above presents us with an easy way to analyse many possible disk-halo 
decompositions. To obtain the core radius and central density of a flattened dark halo, one uses the 
equations ( |B18|) and (B19) and multiplies these values by C{q) (equation [[A^ ) and TLi^q) (equation [|A9t| ) 



respectively. 

B.4. The Edge of the Optical Disk 

Combining the relations for the M/C -ratio, the DM-halo core radius and central density it follows 
that : 



Pstars (^) 

Ph{r) 



27.7 

V 10 



(^c(7) 



2 , 2 




(B20) 



where r = R/hff . Thus the ratio of stellar to DM midplane densities is independent of the amplitude 
of the observed rotation curve and depends only weakly on the shape of the rotation curve (via f3x), 
the truncation radius and thickness of the stellar disk (via fm and fx and hji/ze)- The dependence of 
Pstarsi^)/ Phi^) upon galactoccntric radius is graphically presented in Figure |l4| for several 7 values (full 
line 7=0.9, dotted line 7 = 0.8, short dashed line 7 = 0.7, long dashed line 7 = 0.6, dash-dotted line 
7 = 0.5). I term the distance at which the disk densities no longer dominates^ the total midplane densities, 
Rdde- The maximum-disk like models have Rdde values ranging from four to seven, corresponding to the 
range where some stellar disks are observed to be truncated (vdKS81a&:b; vdKS82a&b). Furthermore 
we notice that the luminous-to-dark matter ratio decreases rapidly with decreasing 7 and q values. 



i.e. become smaller than a few times the DM density. 
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Recently Barteldrees &: Dettmar (1993) have found that some spirals have edges in their light distri- 
bution as close in as two radial scale-lengths. Only for substantially rising rotation curves (say (3 > 1.3) 
and/or small 7 values does Rdde become as small as two, suggesting that these systems have rising rotation 
curves and/or are DM dominated. 

C. A self-consistent solution for the mult i- component method 

In this Appendix, from the equation of hydrostatic equilibrium (|l^) , I derive a self-consistent analytical 
solution for the vertical distribution of the gas for the multi-component method proposed by Bahcall 
(1984). The solution (equation [P28|| ) is in the form of an integral equation for z{pgas), which may be 
inverted numerically to yield Pgas{z)- This equation can be used, in an iterative manner, to determine the 
local dark matter density. The model galaxies presented in the Figures |, H and m were calculated using 



this analytic method. As it turns out, an analytic solution can be found if one makes the approximation 
that the DM-halo density is constant with height above the plane in the region where the gas is co-spatial 
with the DM-halo. For DM-halos which are not too flattened this is not a bad approximation for most 
galactocentric radii. For the gas as well as for the stars the following equation holds : 

crli ^^'^^^ = -AtiG r {pgas{z') + Pstars{z') + {phalo + Prot)} dz' . (C21) 

a z Jo 

The density and dispersion {pi and c^J can represent any isothermal component i, stars or gas. Since 
the right hand side of (C21) is the same for stars and gas, it follows that the vertical distributions are 



powers of each other. The density distribution has to be normalized such that the integral equals the 
surface density. So I write : 



Pgasiz) = Cg [pstarsT^ = Pgas{0) (^£^^1) 

XPstarsK'J)/ 



Psg 

(C22) 



with 



Psg = -A^ and Cg such that Sgas = / Pgas{z') dz . 

^gas J-oo 

Similar relations can be defined for p stars- Differentiating ( |C21D with respect to z yields : 
dj^ In p asiz) 

Agas = Pgasiz) + Pstarsiz) + iphalo + Prot) , (C23) 

with 

A,as = (C24) 
Solving for the gaseous component after defining aiz) = In pgasiz) I rewrite ( p23| ) : 
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A 



gas 



d? a{z) 



e" + Cs (e")P^» + {phalo + Pr. 



ot 



(C25) 



Multiplying ( |C25| ) by 2^ and using the equality : ^ (^)^ = 2 



^ this equation reduces to : 



d f da\' 



2^ {e'- + Cs{e'^Y^^ + {phaio + Prot)} 
dz 



Integrating the RHS with respect to a yields : 



For any physical system, the restoring force vanishes at the midplane (i.e. the RHS of [ C21[ equals zero) 
0. Thus the integration constant C is given by : 



so that ^ 



2 = 



C = -Pgas(O) { {l+p',^gPsg) + {p'h,g + p'r,g) Pgasi^) } , 

where the primed variables are defined as follows : 



(C26) 



s,g 



PstarsjO) . / _ Phalo{0) . / _ ProtjO) 
Ph,g ~ ^ j-Q-j ' Pr,g 



PgasiO) 



PgasiO) 



(C27) 



Thus 



A. 



~ \ ~ \ Pgas + — {Pgasf"" + (Phalo + Prot) ^^Pgas + C 



+0.5 



gas 



Pgs 



da 
dz 

Rearranging terms and integrating yields : 



'Agas [P^-^^'^ dpga 



"2 J PgasiO) Pgas 



c 

Pgas H {PgasY"" + {Phalo + Prot) In /Ogas + C 

\ Pas J 



-0.5 
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Z{pgas{z)) . 



Normalizing the gas distribution, i.e. defining y 
(|C22|) and (|C24|) I find : 



PgasiO) ■ 



inserting the integration constant C and using 



^gas 



2 



\ 8^G Pgas{Q) 



Y y + p's,g Psg) - (y + P'sg Psg V^^') " ip'h,g + Pr,g) Hv) , 



(C28) 



Notice that the term {^p'g^g Vsg) equals the ratio of stellar and gaseous kinetic energy in the midplane. 
For vanishing stellar, halo, and rotational densities (i.e. p'^ g = p'^ ^ = p'^g = 0.0 so that the gas is 
fully self-gravitating), it can be shown that ( p2g| ) then reduces to the familiar self-gravitating solution, 
Pgas{z) oc sedi^ {z / zq) , indeed. 

Investigating equation (C28), we see that there is no solution possible when the sum of the rotation 



and DM-halo density is negative (i.e. for a steeply rising rotation curve) because at small fractional 
densities, Y, the logarithmic term goes to — oo and hence dominates, leading to a negative term inside 
the square root. Bottema et al. (1987), applying the local approach, find that the rotation term can 
contribute significantly in the steeply rising parts of galactic rotation curves. Obviously, since galaxies 
do not have holes in those regions where the rotation curve rises, at least one of the assumptions must 
break down. Most likely, deviations from cylindrical symmetry become important and the equation of 
hydrostatic equilibrium is no longer a good approximation to the vertical Jeans equation. 

An iterative procedure is followed to find the vertical density distribution of the gas. Starting from 
initial guesses for stars and gas, a solution of ( |C28 ) is calculated, which provides a new (normalized) 



z-distributions for the first component. Note that this equation can be inverted to calculate the vertical 
distribution of the stars from the z-distribution of the gas merely by interchanging the g and s indices. 
In order to eliminate numerical instabilities, I chose the calculated component (i.e. the LHS of ]C28|] ) 
to be the component with the smallest velocity dispersion. The normalized distribution, Zg^s (Y) is then 
inverted to yield pgas{z) which is then scaled so that its z-integral yields the appropriate surface density. 
The stellar density distribution is then calculated by raising the gaseous distribution to the power pgg 
(i.e. applying the equivalent of equation |C22(| ) followed by the surface density normalization. From the 
new midplane values the next solution of (|C28| ) is calculated, ... until consecutive changes become small. 
For a given stellar velocity dispersion, the stellar vertical density distribution changes in each iteration 
step. In order to keep the thickness of the stellar disk fixed, as suggested by observation ( vdKS81a&:b; 
vdKS82a&b; BD93), I change the stellar velocity dispersion after each iteration step accordingly. This 
procedure can be easily expanded to include "N" isothermal components. 

The local dark matter density can be determined in an iterative manner as well. An upper limit to the 
local dark matter density can be found by assuming that all other components are massless. Using the 
Gaussian-equivalent method (equation ||l^), the measured width and velocity dispersion of the gas layer 
the upper limit of the dark matter density is calculated. A lower limit to the local dark matter density 
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is of course zero. For a values of the dark matter density halfway between the lower and upper limit, the 
model gas-layer thickness is calculated and compared with the observed value. If the calculated thickness 
corresponding to this trial DM-density value is larger than the observed value, then the lower bound is 
increased to the trial value. This procedure is repeated until convergence is reached. Model calculations 
show that beyond the optical disk, it is possible to accurately recover the input (into eqn. [ C28[ ) DM 
densities from simulated thickness measurements. 

D. Four Special Cases 

Here I describe the gaseous distribution in the vertical direction for four interesting limiting cases with 
very different flaring behavior. In the first three of those, the rotation curve is assumed to be constant, 
while the fourth case represents gas in a region where the rotation curve is falling. In all but the self- 
gravitating case, the gas is treated as a massless testparticle component. The results are listed below, 
where all widths, Wgas{^ FWHM/2.35, as defined by eqn. |5|), are expressed in kpc. 



Fully self — gravitating gas : 



6.74 



10 



(D29) 



Gas in a stellar disk 



0.08, 



^0,400 



X o"io exp 



R 
2h^ 



(D30) 



Gas in a (flattened) dark halo mass distribution : 




2.436 



1.436 + q 



<7lO 



V{CC; l)halo,100 



lRl,, + {R^o/C{q)f 



(D31) 



Gas in Keplerian falloff region : 

o"io 



= 0.484 



totAl 



-"lO 



(D32) 



where aio is the gaseous velocity dispersion in units of 10 km/s, Ei(i?) the gaseous surface density in 
units of solar masses per pc^, 2:0,400 the scale height of the stars in units of 400 pc, /iioo(O) the central 
surface brightness in units of 100 Lq/pc^, Mj L2 is the number of solar masses per solar luminosity, /ir 
is the radial scale length of the (exponential) stellar disk, y(oo; l)/ia/o,ioo the asymptotic rotation velocity 
of the round dark halo in units of 100 km/s, i2c,io the core radius of the dark halo in units of 10 kpc, i?io, 
the galactocentric radius in 10 kpc, C((;), is the core radius correction factor arising from the flatness, g, 
of the DM halo (is of order unity, see equation | A8|), and Mtot,ii the total mass of the galaxy in units of 
10^1 Mq. 
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Notice that the thickness of the gas layer in a potential which is completely dominated by the rotation 
density is of the same order of magnitude as a gas layer in a DM-halo potential only. Such situations 
may arise for galaxies where the rotation curve falls in a close to Keplerian manner. 

If one really wants to use the local approach (e.g. for quick estimation purposes) instead of the global 
approach, one might as well use the approximation to equation ( |C28 ) I found to be accurate to about 
10% over a wide range of parameters. Wg^^^°^ can be found from : 



jx/o-pprox 
, rr gas 



1 



where 

\ 2 



no— self 



gas 



W, 



self 



gas 



yt/stars 
^ ''gas 



+ 



+ 



W-no—sel f 
.gas , 



I Whr 

\ VVhalo 
\ gas 



(D33) 



(D34) 



with 

Wg 

Whr 

and 

'^halo 



(1.15 Sgas + 1-55 {Yistars + ^halo + '^rot)) 
{'^gas + ^stars + ^halo + '^rot) 

(1.20 T,stars + 1-05 i^halo + '^rot)) 
i^stars + '^halo + ^rot) 



Phalo X W, 



no— self 
gas 



and T,rot = Prot x 



no— self 
gas 



where the relations for the halo and rotation surface densities are evaluated (eqn. [D34]) using Whr = 1-0. 
The width s, W^ g^g, are the widths as found for the special cases ( according to the equations | D29 |, | D3C ], 
[ D31 | and | D32| ]). Inside the stellar disk, the width of the gas layer is mostly determined by the gas-inside- 



stellar-disk term, Wg^f"^. Beyond, the halo-term, W^j°, dominates. The self-gravity of the gas perturbs 
the dominant term slightly, but with a different numerical coefficient for the two regimes because of the 
different form of the un-perturbed z-distribution functions (sech^^=9(2;/2:o) versus Gaussian respectively). 

I have tested this approximation for several cases. The first case corresponds to BE89's mass distri- 
bution for NGC 3198. The second case corresponds to a 7=0.6-disk in a flattened {q = 0.1) dark halo, 
thereby minimizing the stellar and maximizing the dark halo contribution. The third case consists of 
BE89's halo and stellar mass model but with the gas surface density multiplied by a factor of ten. In 
the last test case I also multiply the halo densities by a factor of ten. In all these cases equation ( p33| ) 
recovers the exact solution to equation ( C28| ) to within 10 %. 
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Fig. 1. — The lower panel shows the rotation 
curves of a family of DM-halo models with vary- 
ing flattening {q = c/a = 1.0,0.7,0.3,0.1), cen- 
tral density {phfl{q), equation [lAStl ) and core ra- 
dius {Rc{q), equation [|A8[| ). The fact that all 
curves fall on top of each other illustrates that 
the equatorial rotation curve alone is not suffi- 
cient to determine the shape of the DM-halo. In 
the top panel I present the ratio of the rotation 
curve of the flattened to the round DM-halo for 
the same q^s as in the bottom panel. Although 
the residuals show systematic behavior, the am- 
plitudes are smaller than the routinely obtained 
observational uncertainties (i.e. BE89) 



Fig. 2. — In this figure I present the q- 
dependence of the DM-halo central density (top 
panel, equation lAE]) and core radius (lower 
panel, equation [|A8 |). With these functional 
forms for Ph,o{<l) and Rdg), the resultant equa- 
torial DM-halo rotation curve is almost indepen- 
dent of q (see Figure I) . 
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Fig. 3. — Here I show the dependence of the DM-halo parameters on the mass-to- hght ratio of the stellar 
disk. As an example I used the parameters valid for NGC 3198 (L(0)=207.1 Lq/pc^, hji = 2.3 kpc, 
Ze = 0.23 kpc). The gaseous surface density distribution and rotation curve (Vo6s(2.3/i/j) = 145.8) were 
taken from BE89, who used Kent's (1987) photometry to calculate the rotation curve due to the stellar 
disk {fm = 0.899 and /lo = 0.578). In the lower right panel the Ad/jO -ratio is presented as a function 
of 7 (defined by eqn. yB15U ). In the upper left and lower left panel I show the DM-halo central density, 
Phfi, and its core radius, Rc respectively. The dark-to-stellar mass ratio is presented in the upper right 
panel. To indicate the sensitivity of the halo parameters on the slope of the rotation curve, I calculated 
them for three values : the appropriate /Js-value for NGC 3198 (0.985) and the it 5% values. 



Fig. 4. — In this figure I present fits to the rotation curve of NGC3198 with different 7 values. In the lower 
panel I present the observed rotation curve (open squared), the rotation curve due to the Hl-disk (crosses) 
and several disk-halo combinations [ full line, gas+(7=0.9-stellar-disk)+DMhalo,A^/£ =2.96, phfl = 
11.2, Rc = 6.26 ; dotted line, gas+(7=0.7-stellar-disk)+DMhalo,7W/£ =1.79,phfl = 92.6, Rc = 2.00 ; 
dashed line, gas+(7=0.5-stellar-disk)+DMhalo,A4/£ =0.92, phfi = 684, Rc = 0.73 ], where the central 
densities are expressed in units of uiMq /pc^ and the core radii in units of kpc. The upper panel shows 
the velocity difference between observations and model. The large 7 = 0.9-value corresponds to the 
"maximum-disk" solution. A better fit is the 7 = 0.5-solution (low M/C -disk) close to the value pro- 
posed by Bottema (1993). Furthermore, since the velocity field of the inner parts of NGC 3198 deviates 
from axisymmetric rotation (Corradi et al. , 1991), the discrepancies between observed and model curves 
in the inner regions seem irrelevant. The only free parameter in the fitting procedure I used (described in 
Appendix B) is 7 (defined by eqn. B15 |). The observed rotation curve (and the errors, i.e. the difference 
between both sides, are indicated by the error bars), the rotation curve due to the gaseous disk, the 
rotation curve due to the stellar disk and the optical parameters were taken from BE89. 



Fig. 5. — For a model galaxy (consistent with the parameters of NGC 3198 as observed by BE89) we 
compare the midplane mass densities inferred from the thickness of the gas layer {pm-p, using the local ap- 
proach via equation |]l4|, where <T(2 gas)=7.5 km/s) with the true stellar midplane densities (p(STARS)iNP, 
left-hand set of panels), the true (stellar+DM-halo) densities (p(STARS + IIALO)inp, middle set of pan- 
els), and the true (stellar+DM-halo+rotation) densities (p(STARS + HALO + ROT)inp, right-hand set 
of panels). The truncation of the stellar disk at 6/i/j is easily recognized for almost all disk-halo models. 
The self-gravity of the gas has not been included in these calculations. The bottom two panels were 
calculated for a round halo, while the top panels corresponds to a fiattened {q = 0.1) DM-halo. The lines 
correspond to models with different 7 values (full line : 7 = 0.9, dotted line : 7 = 0.8, short dashed line : 
7 = 0.7 and the long dashed line : 7 = 0.6). For most models, the ratio of derived-to-input densities (ver- 
tical axis) deviates significantly and systematically from unity so that the Gaussian-equivalent method 
can not be used to reliably determine the midplane volume densities (see also § ^.lj ) 
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Fig. 6. — In this collage of model galaxies, each plot represents the vertical force, Kz{z\^,q), as a 
function of height above the plane {z). The top plots correspond to a 7=0.9-disk, the lower two plots to 
a 7=0.6-disk. The left-hand side plots were calculated for a round {q = 1.0) halos; the two right most 
plots for a flat {q = 0.1) halos. The different panels inside each plot correspond to different radii (in 
units of /i/j). All panels have been calculated for a galaxy model consistent with the observed properties 
of NGC 3198 (BE89). The open squares represent the vertical force due to a double exponential stellar 
disk. The dashed line is the vertical force valid for the local approximation. The full line is the sum of the 
vertical force from the stellar disk and DM-halo. We notice that the maximum-disk plus round DM-halo 
combination is the only region in 7, g-space where the stellar disk is dominant over the entire extent of 
the stellar disk { R <^ ). For this combination of 7 and the stellar disk contributes significantly (~ 
20 % at 11 hfi) to the vertical force, even beyond the truncation of the stellar disk (at 6/i/j). Furthermore 
notice, that the shapes of the vertical forces generated by a very flattened DM-halo are very similar to 
the shapes of the vertical force in a round DM-halo, but that the amplitudes are about twice larger. 



Fig. 7. — This figure presents the vertical force in the region beyond the optical disk. Again, the galaxy 
model used, stars+DM-halo+gas, is consistent with the observed properties of NGC 3198 (BE89). The 
full line is the sum of the vertical force due to the stellar disk, the DM-halo and the gas layer. The dotted 
lines represent the Kz{z) due to the flattened DM-halo only (g = 1.0 and 0.3 for the Figures]^ (left-hand 
plot) and 0b (plot on the right-hand side) respectively). The vertical force arising due to the (truncated) 
stellar disk is represented by the short dashed line. The long dashed line represents the self-gravity of the 
gas, which' contribution to the total vertical force can be rather important (30-40%) beyond the optical 
disk. 
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Fig. 8. — In this figure I illustrate the depen- 
dence of the thickness of the gas-layer, Wgas{^ 
FWHM/2.35, as defined by eqn. §), upon 
it's self-gravity. I used the multi-component ap- 
proximation discussed in Appendix C to calcu- 
late NGC 3198-like mass models with varying 
amounts of gas. For the mass models with the 
differently scaled gaseous surface densities I per- 
formed a separate disk-halo decomposition, all 
yielding similarly acceptable fits to the observed 
rotation curve. The left-hand set of panels shows 
(for q = 0.1, q = 0.5, and q = I from top to bot- 
tom) the flaring of the gas-layer for several gas- 
surface-density scaling factors. The scale factor, 

S, is defined such that T,gas,Model = sT, gas, observed 

; full line , s = 0.01 ; dotted line , s = 1 ; 
short dashed line , s = 5 ; long dashed line , 
s = 10. In the right-hand set of panels I plot the 
thickness ratio of scaled surface density to the 
zero surface density case. From these plots it is 
clear that the effects of the self-gravity are mod- 
erate, but vary systematically with radius for the 
observed gaseous surface densities. For the up- 
wardly scaled surface densities I find almost no 
dependence of the gas-layer width on the halo 
flattening. 

Fig. 9. — In the bottom panel I compare the 
thickness of the gas layer {Wgasi^ FWHM/2.35, 
as defined by eqn. ||5|)) as calculated in the 
global approach (the lines with the filled circles) 
with the thickness as calculated from the multi- 
component approach (the lines). For both ap- 
proaches a fixed sech^(2;/(22;e)) stellar distribu- 
tion, truncated at R/hR=6, was used. Three dif- 
ferent DM-halo flattenings are shown (from top 
to bottom : q = 1.0, full line; q = 0.3, dashed 
line; q = 0.1, dotted line). In the top panel I 
plot the percentage difference between the two 
approaches. The difference between the global 
and multi-component method mainly arrises due 
to small scale structure in the observed rotation 
curve, which has been used to calculate prot- 



Fig. 10. — For a given rotation curve (BE89), I show how the thickness of the gas layer (Wgas(~ 
FWHM/2.35, as defined by eqn. |5|)) depends on the flattening of the DM-halo {q = c/a). The different 
line-types/symbols correspond to different halo flattenings q ; q = 1.0 (fuh hne), q = 0.8 (dotted line), 
q = 0.6 (short dashed hne), q = 0.3 (long dashed hne)), and q = 0.1 (fihed circles). The model curves 
presented were calculated for a gaseous velocity dispersion of 7.5 km/s. The left-hand side figure was 
calculated using a double exponential stellar disk, which is replaced by a truncated stellar disk for the 
middle panel. For the right-hand side plot I included the truncation of the stellar disk as well as the 
self-gravity of the gas. For these particular mass models the gaseous self-gravity influences the model 
widths more strongly than the truncation of the stellar disk does. 



Fig. 11. — This figure shows how the thickness of the gas layer (Wgas{^ FWHM/2.35, as defined by 
eqn. |^)) depends on the mass-to-light ratio of the stellar disk for several galaxy models (truncated- 
stellar-disk-|-roitn(f-DM-halo) consistent with the observed properties of NGC 3198 (BE89). The different 
line types correspond to different mass-to-light ratios of the stellar disk ; {A4/C =2.98, 7 = 0.9) full line, 
{M/C =2.35, 7 = 0.8) dotted line, {M/C =1.80, 7 = 0.7) short dashed line, {M/C =1.32, 7 = 0.6) long 
dashed line. Inside the optical disk, the extreme thinness of the HI disk (in absolute measures) will, 
quite likely, limit its usefulness for the determination of the /£ -ratio of the stellar disk to the nearest 
galaxies. Beyond the optical disk, the thickness of the gas layer becomes virtually independent of the 
mass-to-light-ratio of the stellar disk, so that a determination of the DM-halo flattening is independent 
of M./C . As for Figure 10 the model curves presented were calculated for a gaseous velocity dispersion 
of 7.5 km/s. 



Fig. 12. — This figure illustrates one of the observable effects of a flaring Hl-layer. The apparent gaseous 
velocity dispersion (i.e. the second moment w.r.t. recession velocity) of two galaxy models is shown here 
(the input dispersion was 7.5 km/s). The upper panel was calculated setting the FWHM of the gas layer 
to 100 pc everywhere, while the lower panel corresponds to a model run where the vertical distribution of 
the HI layer was calculated in a self-consistent manner (cf. the procedure outlined in Appendix C for a 
model with q = 1). A flaring disk imprints a clearly different pattern on the apparent velocity dispersion 
map. Since this is a large-scale pattern, its detection should be relatively easy. The model calculations 
were convolved to a 15"xl5" beam, whereas the channels are square and 2.56 km/s wide. For both models 
I used the rotation curve and HI surface density distribution of NGC 3198 (BE89). 



Fig. 13. — This figure shows that the effects of a flaring HI layer are also visible in individual channel 
maps. The top panel shows a channel map {\Vsys — Vchan\= 104 km/s, smoothed down to 30"x30") for 



the thin hydrogen disk described in Figure 12. In the middle panel I show the same channel but now for 
the flaring-disk model. The "butterfly-wing" for the thin hydrogen layer model (upper panel) is clearly 
narrower and more peaked than its flaring equivalent (middle panel). Because the difference between the 
two models (lower panel) is significant over a large range in radii, it is expected that local irregularities 
(in the velocity field and/or the surface density distribution) will not be a limitation in determining the 
thickness of the hydrogen layer of inclined galaxies. 
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Fig. 14. — In this figure I present the radial dependence of the ratio of stellar to dark matter midplane 
density for several 7 values (full line 7 = 0.9, dotted line 7 = 0.8, short dashed line 7 = 0.7, long dashed 
line 7 = 0.6, dash-dotted line 7 = 0.5). The left most figure corresponds to a round DM-halo {q=1.0), the 
right most to a very flattened DM-halo {q=0.1) and the figure in the middle was calculated for a DM-halo 
of intermediate flattening (g=0.4). This ratio is calculated using the disk-halo decomposition for a flat 
rotation curve presented in Appendix B (| B20| ). It can be clearly seen that in maximum-disk like models 
the stellar density dominates the midplane mass density out to about four radial scale-lengths, beyond 
which the contribution of the stellar density decreases sharply. For the low 7- models, the stellar density 
never really dominates. 
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